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Abstract 

We consider a dilute solution of infinitely rigid rods near a curved, perfectly 

repulsive surface and study the contribution of the rod depletion layer to 

the bending elastic constants of membranes. We find that a spontaneous 

curvature state can be induced by exposure of both sides of the membrane 

to a rod solution. A similar result applies for rigid disks with a diameter 

equal to the rod's length. We also study the confinement of rods in spherical 

and cylindrical repulsive shells. This helps elucidate a recent discussion on 

curvature effects in confined quantum mechanical and polymer systems. 
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The elastic properties of fluid membranes are among the ruling factors in the physics 
of many biological and surfactant systems such as cell membranes, vesicles, lyotropic liq- 
uid crystals or microemulsions ]IJ. When macromolecular species are also present in such 
systems, the interactions between the membranes and the macromolecules change the bare 
elastic constants. Recent studies have explored contributions from adsorption de- 
pletion [[§-[(| or end-grafting °f flexible polymers. In all of these cases, exposure 
of both sides of the membrane to a solution results in a modification of the bending and 
splay modulii K and K, the sign of the contribution depending on the possibility of surface- 
bulk equilibration. For instance, depleted or adsorbed polymers reduce K and increase K, 
whereas the reverse holds for end-grafted polymers. 

Interest in solutions of rigid macromolecules can be traced to the extraction of tobacco 
mosaic virus (TMV) [|Kj and to the subsequent observation of the nematic phase in TMV 
solutions PHI , later explained by the seminal theory of Onsager fll2f . Rod shaped particles 



in the colloidal range have since been studied in a large variety of mineral and organic sys- 
tems [|nj , they are also present in the biological realm where examples range from TMV-like 
virus to fibrils of amyloid /3-protein, the molecular agent at the origin of the Alzheimer dis- 
ease. The surface interactions of rigid macromolecules were first studied by Asakura and Oo- 
sawa (T^j who showed that the steric depletion of the molecules at a flat surface increases the 



interfacial energy of the system, implying that two surfaces brought to separations smaller 
than a rod length will experience an attractive force. The theory of depletion interactions 
for rods have since been extended to include ordering effects of the bulk phase |16|] or effects 
of rod-rod excluded volume interactions | 17| ; but it has been mostly |18|] devoted to flat ge- 



ometries. In this letter we extend the original Asakura and Oosawa theory to include effects 
of surface curvature and determine the contribution of the rod depletion layer to the elastic 
constants of the membranes. An important reason for considering this case can be better 
understood by comparing the order of magnitudes of changes in K and K caused by depleted 
layers of rods with those caused by spherical particles. In a solution of colloidal spheres of 
radius r and particle number density pb, the typical scale of the energy density is ksT pb- 



Corrections to the interfacial energy are thus of the order A7 ~ ksTpbro- On other hand, 
interfacial tensions in most liquids are of the order of 70 — ksT/a 2 where a is a microscopic 
size. For instance for a ~ O.lnm, 70 is of the order of tens of mN/m. The corrections due 
to the depletion of spherical particles are thus a factor (a/ro) 2 lower than typical interfacial 
tensions even at order unity volume fractions <fi = pb(4n /3)tq. Curvature corrections to the 
interfacial energy can generally be written as A7 ~ fcsTp6r (l + C\r /R + C^r^/i? 2 ), with 
R the radius of curvature and C\,C2 two numerical constants. This leads to modifications 
of the bare elastic constants of the order of AK ~ ksTpbrQ-. at the upper concentration 
limit pb ~ l/ r 0) these corrections are only of order ksT, a value at the lower end of the 
range |19| |, 1 — 20k bT, of most bare elastic constants. In the case of a rod solution with rod 
number density pb, the upper concentration limit of the isotropic solution is the Onsager 
concentration pi = 4.2 /(L 2 D) ||12|| , with L the length of the rod and D its diameter. The 
contributions to the interfacial tension are now of the order A7 ~ ksTpbL, but even for 
rod diameters of order of the microscopic length a the contribution to the interfacial tension 
of rod solutions at the Onsager concentration is still a factor {a/L) smaller than typical 
interfacial tension values. However, modifications of the elastic constants are here of order 
AK ~ ksTp^L 3 , a factor {L/D) larger than k^T . Therefore, even rather rigid phospholipid 
membranes with elastic constants as large as 20k bT may have their rigidities modified sub- 
stantially at low rod concentrations where rod-rod interactions are negligible. We remark 
also that the steric surface interactions are identical for rigid rods and rigid disks if one takes 
the rod's length and the disk's diameter to be the same; however the disk system does not 
enjoy the above-mentioned L/D enhancement. 

We consider an ideal gas of rods of length L in the presence of flat and curved surfaces 
which repel the rods. We parameterize the possible rod configurations by the center of mass 
coordinates, r, and the two angles specifying in which direction the rod points, u = (#,</>). 
For sake of simplicity the thickness of the rods is taken to be zero. The relevant potential 
describing the thermodynamics of the system can be written as 
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Q [p(r, uj)} = J df J duj p(f, u) [log(u p(r , u)/e) — 

(fi b -U ext (f,uj))] (1) 

where v is some normalization volume, [i b is the solution chemical potential, and U ext is 
the hard wall interaction potential that is either infinite or zero, depending on whether 
the configuration of the rod is allowed by the 'hard wall' requirement or not. Functional 
minimization of Q with respect to p(f, uj) gives the equilibrium density profile: 
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p(r, uj) = — e - u ^^) = £L e -u ext (f,u) (2) 

V 47T 

The local, position dependent, number density of particles p(r) is simply the sum over all 

allowed angular configurations of p(f,uj). The excess surface energy is then calculated from 

fl [p(f, u)} — Q [p(z — > oo, uj)} 

S = 
dz(p b -p(z))J(z,R) (3) 
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where S is the surface area, z the perpendicular distance from the surface and J(z, R) a phase 
space factor, which depends on the geometry. For flat surfaces J(z, R) = 1, for cylinders 
of radius R, J(z,R) = 1 + z/R, and for spheres J(z,R) = (1 + z/R) 2 . From equation (^) 
it is clear that differences between the excess energy of a flat and a curved surface depend 
on two factors. The first is the configurational part measured by the differences in the rod 
density profiles, the second one is associated with the space available to the center of mass 
in the neighborhood of the surface, and is measured by J(z, R). In the case of hard sphere 
solutions, where there is no coupling between configuration and curvature, only J(z, R) is 
responsible for energy differences; (p(z) = pbQ^z — L/2), independent of geometry). It is 
easy to show that in that case one has A7 = kBTpbr (l + r /(2R)) for cylindrical surfaces 
and A7 = k B Tp b r (l+r / R+r^/ (3R 2 )) for spherical ones. Following the usual procedure U 
one then finds the corrections to the elastic constants of a membrane exposed to a colloidal 
solution: AK = and AK = 2/3 k B T p 6 rg. 

Figure 1 shows the angular geometrical constraints for a rod close to a flat or curved 
surface. When the center of mass is at a perpendicular distance larger than L/2 from the 



surfaces, the direction of the rod is not constrained. For distances < z < L/2 only a 
fraction of the solid angle is available to the rod. For rods outside a curved geometry a 
second distance z° must be defined: for z° c < z < L/2 the rod touches the surface with its 
end, whereas for distances smaller than z° the rod touches the surface somewhere along its 
length. Inside a curved geometry yet a different length z l c must be defined, below which 
the rod is geometrically not allowed. Performing the angular integrations one gets the 
density profiles plotted in figure 2. As expected rods outside a spherical curved surface 
have a larger density than rods close to a flat surface. Conversely, rods inside a spherical 
surface have a smaller density than rods close to the flat surface. The density profile of 
rods close to cylindrically bent surfaces is identical to the spherical profile when the rods 
are oriented perpendicularly to the cylinder axis, and to the flat profile when the rods are 
oriented along the cylinder's axis. Intermediate angles interpolate smoothly between the 
two limits. A notable feature of the density profiles is that they are non-analytic functions 
of the curvature. For instance, the inner part (z < z°) of the external profile depends on 
the square root of the curvature, and the internal profile has a cutoff at z l c . This implies 
that even for large radii of curvature one cannot simply obtain the internal profiles from the 
external one by a R — > —R transformation. Integration of the profiles lead to the following 
contributions to the interfacial energies: 

A7 out = k b Tp b ^ (4) 

for rods outside spheres and cylinders, or close to flat surfaces; 

A 7in = k b Tp b j ^1 - a^j (5) 

for rods inside spheres (a = 1/12) and cylinders (a = 1/32). Results for the outside configu- 
ration are exact. For the inside configuration results are exact for spheres, and perturbative 
to second order in 1/R, for cylinders. Equations (|]) and (|^) bear a few interesting conse- 
quences. For instance equation (|j) implies that the excess free energy of a convex volume 
immersed in a rod solution does not depend on how the surface of the object is curved, but 



only on its area. Also, equation (||) indicates that flexible membranes which expose one 
surface to a rod solution will spontaneously bend towards the solution. Note however that 
there is no spontaneous curvature in the usual sense of a contribution to the energy linear 
in 1/R. When the membrane is immersed in the solution (exposure of both sides of the 
membrane) the total surface energy excess is given by the sum of the two contributions A7 in 
and A7 out . Curving the surface still decreases the total free energy, but in this symmetric 
situation there is no preferred side toward which the membrane should bend. Immersion 
will therefore spontaneously break the symmetry of the system. To the extent that the non- 
analytic character of the total free-energy can be neglected (there is a finite discontinuity 
in the second derivative at 1/R = 0), the following contributions to the membrane elastic 
constants can be extracted 

AK = -k B Tp b — = -k B T^ — -, (6) 
^ 64 B pt 15.2 D" y ' 

AK = k B Tp b — = k B T— — - — — . (7) 
^ 96 p* b 223D v ; 

One finds thus a reduction in K and an increase in K - this favors the formation of saddle 
structures (for instance, cubic phases with periodic minimal surfaces). This is similar to 
other depletion and equilibrium adsorption problems, but as explained above the amplitude 
of the contributions in this system is very large (of order of L/D times larger than k b T, a 
value to be compared for instance to k b ThaN for adsorbed polymers with polymerization 
index N). 

In a confined geometry, e.g. in ordered stacks of membranes (lyotropic smectics), in 
multilamellar vesicles or between any other two surfaces separated by some distance d the 
rod configurations are restricted by two constraints. If the distance d is larger than the rod 
length the confinement effect can be calculated from the above results simply by summing 
the effects of two non-overlapping depletion layers. If the surfaces are at a distance smaller 
than the rod's length the two depletion layers overlap, and one needs to compute two- 
surface effects specifically. The overlap of the rod depletion layers is known to lead to 
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an attractive force between the two surfaces, but recent calculations that include rod-rod 
excluded volume interactions in flat geometries show that a repulsive force can set in at 



distances of the order of the rod length ||17||. Our method can in principle be applied to 



the evaluation of the force between two convex spherical or cylindrical surfaces, but we 
postpone presentation of those results to a forthcoming paper and concentrate here on 
the simpler case of confinement in spherical or cylindrical hard shells. This problem is 
related to a recent discussion of the curvature effects in confined quantum mechanical or 
macromolecular systems. Indeed, if a quantum mechanical particle is confined in a one or 
two dimensional manifold embedded in three dimensions, its movement along that manifold 
will be described by the usual Schrodinger equation with an additional curvature dependent 
potential. For instance a quantum particle moving in a confinement tube will always be 
attracted by curved regions of that tube. A quantum particle in a two dimensional manifold 
will experience an extra potential pOfl , U = — ^) > with Ri and R 2 the two 



principal radii of curvature. Particles will be attracted to curved regions of the surface 
except when it is spherically bent (Ri = R 2 ). It has also been shown that such curvature 
effects propagate into the classical world of statistical mechanics and that the quantum 
mechanical results hold also for flexible polymers confined into spherically or cylindrically 
bent gaps . For such a system the difference in the free energy of confined states in curved 
shells as compared to flat ones is given by U = —k B Tp b R^^ (-^ — -^j , where R 2 G = N a 2 , a 
being the monomer size. Apart from the typical energy scale, these potentials are identical. 
This identity is due to similarities in the equations governing the behavior of both systems, 
and in fact holds more generally for other systems containing Laplacians in their equations 
of motion. In this context it is also interesting to investigate the confinement effect for 
other classical systems, whose configurations are not determined by Laplacian equations. 
For confined colloidal particles it is convenient to study these effects by considering open 
curved gaps, i.e. shells that allow the concentration of their solutions to equilibrate with 
an external reservoir. Attractive curvature effects under these conditions will show up as 
an increase of the average concentration in the curved gap as compared to non curved gaps 



of identical thickness. For illustration we compute first to second order in 1/R the average 
concentration for a confined solution of dilute hard spheres: 

Pflat = Ml - 2 t) 

p sph = Ml-2^)(l-^f) (8) 

p cyl = mi - W 

A spherical curvature is therefore repulsive to the particles, whereas the cylindrical bent 
shells behave as flat gaps. For the confined rod-solution we report the value of the average 
concentration for the flat, spherical and cylindrical geometries when the surface separation 
is equal to the rod length (d = L) 



Pflat(d 


= L) = 


pb 

2 




Psph(d 


= L) = 


Pb 
2 


v 1 16 R 2 > 


Pcyl(d 


= L) = 


Pb j 
2 


^ 128 R 2 



At this separation the spherical shell is repulsive for the rod particles, and the cylindrical 
shell is attractive. For lower distances both shells become repulsive. Results for several 
different systems are presented in Table I. Despite quantitative differences, it is apparent that 
cylindrical shells are always more attractive than spherical shells. Comparison between the 
hard spheres case and the rod and the polymer case also indicate that the coupling between 
curvature and configuration reinforces the primary effect of volume exclusion. Note however 
that all the systems quoted in the Table have a repulsive interaction with the confining 
walls. Further insight into this problem and a test for the generality of the conclusions 
drawn here will require the study of systems where the particles are attracted to or bound 
to the confining walls. Work under progress includes for instance end-grafted polymers and 
adsorbed rod-like particles. 
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FIGURES 




FIG. 1. Configurational space of a rod near a flat surface, outside a curved surface, and inside 
of a curved surface. The rods are represented at the angle where they touch the surface. The 
complete allowed angular space is obtained by rotation of the figure around the z-axis. 
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FIG. 2. Profiles for the the rod number density p(z), in p^ units, close to: a) a flat surface; b) 
outside a spherical surface; c) inside a spherical surface. In the particular case shown here we have 
L = 1, R = 2, z° = 0.061 and z\ = 0.064. 
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FIG. 3. Table of the potentials induced by curvature in several systems 
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